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We study numerically the behavior of continuous-time quantum walks over networks which are topologically 
equivalent to square lattices. On short time scales, when placing the initial excitation at a corner of the network, 
we observe a fast, directed transport through the network to the opposite corner. This transport is not ballistic 
in nature, but rather produced by quantum mechanical interference. In the long time limit, certain walks show 
an asymmetric limiting probability distribution; this feature depends on the starting site and, remarkably, on the 
precise size of the network. The limiting probability distributions show patterns which are correlated with the 
initial condition. This might have consequences for the application of continuous time quantum walk algorithms. 

PACS numbers: 05.60.Gg,05.40.-a,03.67.-a 



I. INTRODUCTION 

The study of quantum mechanical extensions of classical 
transport processes has witnessed a considerable growth over 
the last decade, with much attention being given to quantum 
information and the application of transport processes in po- 
tential quantum computers 1 1]. This led to the development 
of several quantum algorithms specifically designed for these 
quantum computers. Among these algorithms are the so- 
called quantum walks, see for instancejyl 01, which are the 
analog of the classical random walks |4, 5], for an overview 
see |6[. 

Such quantum walks have been studied in different sit- 
uations. For instance, it was shown that continuous-time 
quantum walks can cause a remarkable speed-up of transport 
through certain graphs |7], a feature which is not universal, as 
discussed in Ref. 1 8], where it was shown that quantum trans- 
port can also become much slower than the classical one, de- 
pending on the initial conditions. Another variant of the quan- 
tum walk, the quantum mulitbaker map, was shown to ex- 
hibit a crossover from classical to quantal behavior with time, 
where the crossover time is given by the inverse of Planck's 
constant lot fToll . 

For simple structures these quantum walks are directly re- 
lated to well known problems in solid state physics. Thus, 
methods used in solid state physics can easily be applied to 
quantum walks. For instance, walks on one-dimensional lat- 
tice with periodic boundary conditions are readily treated by 
a Bloch ansatz 

Here, we study quantum walks on two-dimensional net- 
works. We focus on structures topologically equivalent to 
square lattices. For these we evaluate the transition probabil- 
ities between different nodes of the finite networks and com- 
pare these to the Bloch solutions. As we are going to show, 
these transition probabilities display odd, unexpected features. 

The paper is organized as follows: In Sec|II|we briefly re- 
view the properties of continuous time quantum walks and the 
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Bloch ansatz. The time dependent transition probabilities be- 
tween the nodes of the network are presented in SecHH] We 
conclude with a summary of results in Sec lIVI 

II. CONTINUOUS TIME QUANTUM WALKS 

In quantum information theory, qubits on a graph are used 
to define the quantum analog of a random walk. There is a 
discrete |4] and a continuous-time |5'| version. However, dis- 
tinct from classical physics, these two are not equivalent to 
each other. Here, we will focus on continuous-time quantum 
walks on networks. 



A. Continuous time quantum walks on graphs 

We consider two-dimensional graphs, topologically equiv- 
alent to finite, square lattices of side length N. In this way the 
nodes of the graph are connected in a very regular manner. In 
general, to every graph there corresponds a discrete Laplace 
operator, sometimes called the adjacency or connectivity ma- 
trix A = (Aij), defined by letting the non-diagonal elements 
Aij equal — 1 if nodes i and j are connected by a bond and 
otherwise. The diagonal element An is given by the number 
of bonds which exit from node i, i.e., An equals the function- 
ality fi of the node. Thus, in our case we have fi = 2 if node 
i is located at a corner of the square, fi = 3 if the node is 
located along an edge (and is not a corner node), and fi = 4 
otherwise. 

Classically, continuous-time random walks (CTRWs) are 
described by the master equation 1 141 H51 

where Pk,j(t) is the conditional probability to find the walker 
at time t at node k when starting at time at node j. We 
assume an unbiased CTRW such that the transmission rate 7 
of all bonds are equal. Then the transfer matrix of the walk, 
T = (Tkj), is related to the adjacency matrix by T = — 7A. 
Formally, Eq.Q is solved by 

Pk,j(t) = (k\e Tt \j). (2) 
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The quantum mechanical extension of a CTRW, the 
continuous-time quantum walk (CTQW), is now defined by 
identifying the Hamiltonian of the system with the (classical) 
transfer operator, H = — T, (3, Q, |8|] . However, there is no 
unique way of defining a CTQW. As mentioned in |2], there is 
some freedom in defining the diagonal elements of the Hamil- 
tonian for a certain graph. A requirement for the CTQW is 
that H define a unitary process, whereas classically the CTRW 
is a probability conserving Markov processes, which requires 
Sfc^fcj = 0- As mentioned in [2], for regular networks, 
where all nodes have the same functionality, different choices 
of the Hamiltonian give rise to the same quantum dynamics. 
Nevertheless, in what follows we directly identify the Hamil- 
tonian with the transfer operator since some of the networks 
we consider are non-regular. 

In the spirit of a localized orbital picture we now introduce 
the states \j) which are localized at the nodes j of the graph, 
requiring in an obvious notation that (k\j) = S k ,j and that 
J2j = 1 where 1 is the identity operator. Under these 

conditions the Schrodinger equation (SE) reads 



.d 



H|j>, 



(3) 



where we set h — 1. Starting at time to from the state 
\j) the system evolves as U(i, to)\j), where U(i, to) = 
exp(— «H(t — to)) is the quantum mechanical time evolution 
operator. Hence, the transition amplitude a k ,j(t) from state 
\j) at time to state |fe) at time t is 

(k\e' mt \j). 



a k,j (t) 
With Eq.Q we find that 
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We note that the main difference between Eq.© and Eq.@ 
is that classically J2 k Vk,j (t) — 1, whereas quantum mechan- 
ically J2k \ a k,j{t)\ 2 = 1 holds. 

For the full solution of Eqs.Q and l|5} all the eigenvalues 
and all the eigenvectors of T = H (or, equivalently, of A) 
are needed. Let A„ denote the nth eigenvalue of A and A the 
corresponding eigenvalue matrix. Furthermore, let Q denote 
the matrix constructed from the orthonormalized eigenvectors 
of A, so that A = QAQ -1 . Now the classical transition 
probability is given by 

(felQe-^Q-^j), 



Pk,j(t) 



(6) 



whereas the quantum mechanical transition probability is 

n k)j (t) = \a k j(t)\ 2 = KfelQe-^Q- 1 ^)! 2 . (7) 



The unitarity of the time evolution operator prevents the 
quantum mechanical transition probability from having a def- 
inite limit when t — » 00. In order to compare the classical 
long time probability with the quantum mechanical one, one 
usually uses the limiting probability (LP), i.e., the long time 
average of irk,,j(t) 11611 : 



Xk,j 



lim — 



dt n k ,j{t). 



(8) 



B. Boundary conditions and the Bloch ansatz 

Given that our two-dimensional structures have side length 
N, they contain N 2 nodes giving rise to N 2 basis states. We 
switch to a pair notation, by which we set \j) — \j x ,j v ), 
where j x and j y are integer labels in the two directions, with 
jx, jy £ [lj N], see FigQ This labeling of the states is not to 
be confused with the labeling of the adjacency matrix. Note 
that capital bold letters denote matrices, while small bold let- 
ters denote the nodes and the states. 
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FIG. 1 : Sketch of a square network arranged as a regular lattice with 
the appropriate numbering of the nodes. Note that the actual geomet- 
rical realization can be much more flexible, see text for details. 

We stop to recall a basic fact, namely that our matrix A 
does not necessarily imply that the structure we consider obey 
any kind of translational symmetry; in fact the symmetry of 
A is topological in nature, a fact well realized both in poly- 
mer physics ifTvtl and also in quantum chemistry lUsll . where 
the A-matrices are fundamental in Hiickel molecular orbital 
calculations [19, 20]. Regular, translational invariant lattices 
are only one possible realization of the network structures we 
are facing here. The methods of solid state physics fL?! fl3ll 
apply, however, even in our very general case. 

The focus in solid state physics is on systems where Born 
- von Karman periodic boundary conditions (PBC) are as- 
sumed. Here, all transmission rates are taken to be equal. 
Hence, time is given in units of the inverse transmission rate 
7 _1 or, equivalently, we set 7 = 1. Now, for an internal site of 
our network (not on an edge or in a corner), the Hamiltonian 
acting on a state \ j) = \j x ,jy) reads 



2 \jx,3y) 
2 \jx,jy) 



\jx + 1, jy) 
\jx,jy + 1) 



\jx ~ l,jy) 
\jx,jy-l).(9) 



PBC extend this equation to all the sites of the network by 
interpreting every integer j x and j y to be taken modulus N. 
With this generalization, the time independent SE 

H|* e ) = E e \^e) (10) 
admits (as is well known) the following Bloch eigenstates 
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N 

N ^ 



e-^')|j>. 



(11) 



as solutions, where 9 ■ j stands for the scalar product with 
9 = {6 X , 6 y ). The usual Bloch relation can be obtained by 
projecting \^e) on the state \ j) such that ^fg(j) = (j\^e) = 
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a-W-V/N, thus * {j x + + 1) = e-^+^^eO'^,^)- 
The PBC restrict the allowed values of 0. In our case (side 
length AO, the PBC require that * fl (AT + l,j y ) = Vg(l,j y ) 
and ^fg(j x , N+ 1) = ^g(j x ,l). It follows that one must have 
9 X = 2nn/N and 9 y = 2lir/N, where n and I are integers 
and n,l £ [1, A 7 ]. It is now a simple matter to verify that the 
also obey (* e |4v) = Ve' and Ee l*e)(*e| = 1. 
Moreover, by inverting Eg. (II 1> one has 

li) = iE ei(e ' i) l^)' < 12 > 

and l£) might be viewed as a Wannier function of the problem 

Qui. 

Furthermore, from Eas.( llO> and dill the energy is obtained 

as 

E e = 4-2cos6 x -2cos6 y = E ea: +E ey , (13) 

with Eg x = 2 — 2 cos 6*2 and Eg y = 2 — 2cos# y . Under 
PBC the two-dimensional eigenvalue problem separates into 
two one-dimensional problems. 

The transition amplitude at time t from state \j) to state |fe) 
is now, using Eq.lQj} twice: 

e,e' 

= J_£y«M e -w-(fc-i) (14) 
e 

In the limit N — > oo, the sums in Eq. J14t may be changed to 
integrals; by making use of Eq.(Q3J we obtain 

7T 

— iAt p 

lim a fc? -(t) - ^—5- / d9 T e- l9 ^- j ^e l2tcose - 

N^oo ,JV 47T 2 J 

— 7T 

7T 

x y dOye^^-^e 12 ^ 03 ^ (15) 

— 7T 

= f fc. -i. -i v e -m j K _. m {2t) Jky _ h {2t) t 

where J n (x) is the Bessel function of the first kind [21]. Thus, 
on a network topologically equivalent to a square lattice with 
PBC the transition amplitude between the nodes j and k is 
given by 

lim TT kd {t) = [J k ^ 3x (2t)J ky _ 3v {2t)] 2 . (16) 

N — >oo 

For systems without PBC the situation is more subtle. If we 
are interested in the behavior of a particular network of finite 
size N x N we have, in general, to rely on the (numerical) 
solution of the eigenvalue problem and on the calculation of 
the corresponding transition probabilities. However, for the 
smallest network of 2 x 2 nodes, which is equivalent to a line 
of 4 nodes with PBC, the analytic result for tt^j (t) was given 
in lHUl [see Eq. (19) there]. There, the Wk j(t) were simple 
trigonometric functions. For larger networks with PBC, the 



TTk, j (t) are given by a multiple sum which follows from the 
Bloch ansatz, see Ea. J14> . Thus, analytic results for large fi- 
nite networks are hard to get. The situation becomes even 
more complex if there are reflecting boundaries. Evidently, 
the expectation is that larger systems will still behave as in the 
PBC case, since for very large systems most of the nodes are 
far from the boundaries. 

The discrete version of the quantum walk on one and two 
dimensional networks was studied in 12211 . Comparable to the 
periodicity found for the CTQW on a line of 4 nodes in filll . 
it was found that in this case also the discrete quantum walks 
show periodic behavior, see also Il23ll . In a related context, 
mixing properties of discrete quantum walks have been stud- 
ied in I I24II25II . Here again, it was shown that a line of 4 nodes 
has a uniform mixing property, i.e. at certain times the tran- 
sition probabilities Ttk,j(t) f° r a U nodes k are equal. This is 
also true for the CTQW on a line with 4 nodes lllll . 



III. PROBABILITY DISTRIBUTIONS 

In the remainder of this paper we calculate for several finite 
networks of various sizes the transition probabilities tt^j and 
their long time average Xk,j- The numerical determination 
of the eigenvalues and of the eigenvectors proceeded using 
the FORTRAN eigensystem subroutine package (EISPACK) 
10. 



A. Eigenvalue spectra 

In order to have a comparison to our calculations, we con- 
trast the numerically obtained eigenvalue spectra for finite net- 
works to the results obtained for PBC from the Bloch ansatz, 
Eq.(ll3i. We note from the start that Eq.J13i limits the latter 
eigenvalues to the interval [0, 8]. Furthermore, the eigenvalues 
and 8 are nondegenerate because they are only obtained for 
n = I = N and for n = I = N/2, respectively. 
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FIG. 2: Eigenvalues X„, arranged in ascending order, for the net- 
works of size N = 9 (black) and N = 15 (grey) obtained for the 
finite network, solid lines, and from the Bloch solutions with PBC, 
dashed lines. 

Figure [2] shows the eigenvalue spectra for two networks 
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with sizes N = 9 and TV = 15 each. Displayed are the results 
for the finite networks obtained by numerical diagonalization 
and for the PBC cases from the Bloch solutions. The numer- 
ical results for PBC agree exactly with the Bloch solutions 
within the precision of our calculations and would be indistin- 
guishable from each other in PigEJ therefore we display only 
the former. We note that the numerically determined spectra 
are also bound by and by 8, although the value 8 is only 
approached in the limit N — > oc. 

Figure[5]displays the findings for a larger network, N = 60, 
again comparing the spectra of the finite network with those 
obtained from Eq. (I13> . We again note that the spectra get to 
be very close (a fact which was to be expected because now N 
is larger than in Fig|2}, but that the PBC case still lies system- 
atically higher. To highlight this fact, we plotted in the inset 
of Fig|5]for each n the difference in the values between the A„ 
obtained with and without PBC. 




n 



FIG. 3: Eigenvalues \ n , arranged in ascending order, for the network 
of size N = 60. The numerical result for the finite network, solid 
line, is compared to the Bloch solution with PBC, dashed line. The 
inset shows the difference between these two results as a function of 

n. 



B. Special initial conditions 




time t 

FIG. 4: Probabilities for a CTQW to be at the initial site at time t on 
networks of size (a) N = 9 and (b) N = 15. The starting sites are 
the middle node m (solid lines) and the corner node c (dashed lines), 
respectively. The results for the finite network, obtained by using 
Eq.0, solid and dashed lines, are compared to the Bloch solutions 
with PBC, obtained via Eq, jl4t . dots, see text for details. Time is 
given in units of the inverse transmission rate 7 -1 . 



given that now the boundaries play an important role. 
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We now turn to the dynamics of the propagation through 
the network and focus on special initial points. In doing so, 
we continue to compare our results for the finite network to 
the Bloch solutions. The computations for the finite network 
without PBC are performed by calculating all eigenvalues X n 
and all eigenvectors of the matrix A using the FORTRAN 
EISPACK routine. The Bloch solutions are obtained by us- 
ing the standard software packages MAPLE 7. 

We begin by choosing the middle node m of the network 
as the starting site. Figure |4] shows the probability of being 
at m after time t obtained via Eq.0 for two network sizes, 
N = 9 and N — 15, respectively. Moreover, we have also 
inserted into Fig0]the results obtained by starting at a corner 
node c. Clearly, on short time scales our numerical results for 
^m,m(t), solid line, agree nicely with the Bloch solutions, 
dots, given in Ea. (ll4> . even for the N — 9 network. However, 
the probability Tr cc (t) of being at the initial corner node c after 
time t, dashed line, differs quite early from the Bloch solution, 



FIG. 5: Comparison of the probabilities 7r m , m (i) for network sizes 
N — 15 and N — 35 with the bulk system where N — > 00. Time is 
given in units of j^ 1 . 

In order to show for finite networks the deviations from 
the Bloch solutions, we plot in Fig|5] for networks of sizes 
N = 15 and TV = 35, the probabilities Tr mtm (t) to be at time 
t at the middle node m, while also starting at m. On short 
time scales, the 7T m>m (t) behave exactly as the Bloch solu- 
tion for the bulk system given in Eq.( ll6t . However, for finite 
networks at longer times there is constructive interference due 
to the reflections at the boundaries, which result in a higher 
probability to find the walker back at m. This interference, 
of course, can only take place after the wave has crossed the 
whole network and returns to m. For the two examples given 
in Fig|5] this happens approximately at the times jt « 6.1 
for TV = 15 and jt « 15.6 for N = 35. As for the one di- 
mensional case studied in fTHl . such deviations are found in 
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general at times around N/2. 
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FIG. 6: Probabilities for a CTQW starting at one corner node c to be 
at the opposite corner node oc at time t on networks of size (a) N = 
9 and (b) N = 15. The results for the finite network, obtained by 
using Eq.Q, dashed lines, are compared to the Bloch solutions with 
PBC, obtained via Eq. <14> . black solid lines. We further compare 
also to the Bloch solution with PBC, dotted lines, where the initial 
and final node is the node k, see Fig|4|and text for details. Time is 
given in units of the inverse transmission rate 7 _1 . 

In Fig|6]we show the transition probabilities n oc ,c{t) to go 
from one corner node c to the opposite corner node oc in time 
t. Again we take N — 9 and N = 15 as network sizes. 
We remark that already after a short period of time there is 
a considerable probability for the CTQW to be at the oppo- 
site node. For example, for the N = 9 network we find that 
ftoc,c(t) ~ 0.16 at jt w 5.4. This means that for jt sa 5.4 
a remaining probability of 0.84 is distributed among the other 
80 nodes, which is roughly an order of magnitude less than 
0.16. 

Not only is there a very high probability to go to the oppo- 
site node, but also is the transport to this node very fast. The 
same holds for the N = 15 network. Here, the first peak of 
^oc.c(t) occurs at about jt w 8.4, for which 7T OC)C (t) w 0.08, 
which is again a relatively high value. We also remark that 
the shortest "chemical" distance between two opposite corner 
nodes on this network is Ax — 28 bonds. Therefore, the (ini- 
tial) "velocity" of the CTQW is Ax/jt « 3.3. 

For a closer examination, we also confront our calculations 
to the Bloch solution, see Fig|6] The results from the Bloch 
ansatz naturally compare only to ours where the middle node 
is the initial node. However, we note that the maxima of the 
probability 7r OC)C (i) on the finite networks without PBC from 
a corner c to the opposite corner oc occur at approximately 
the same positions as the ones for Hk,k{t) from the Bloch so- 
lution for CTQWs with PBC from any node k to the same 
node in the same time t. This is quite remarkable because the 
distances traveled by the CTQWs are different in both cases. 
For the finite network of size N — 15 without PBC, there are 



28 bonds from c to oc, whereas for a network with PBC there 
are always 15 bonds from k to a nearest site corresponding to 
the same k; for the N — 9 network, the corresponding dis- 
tances are 16 bonds and 9 bonds, respectively. This implies 
that in this particular situation the initial "velocity" for the fi- 
nite network is higher than the one for the network with PBC. 



7Tfc,c(t) 7Tfc,c(*) 




FIG. 7: Snapshots of the probabilities TTk, c to be at time t at node 
k when starting at node c = (1, 1). Here the network is of size 
iV = 15 and time is given in units of 7 _1 . 

In FigQwe display snapshots in time of the transition prob- 
abilities TTk, c (t) to go from the corner node c = (1, 1) to the 
other nodes. We see again that on short time scales the trans- 
port is very fast and, furthermore, that the main fraction of 
the probability stays on the diagonal. One might think that 
the effect is ballistic, but this is not the case: neighboring sites 
along the diagonal are not directly related via A (one may also 
remember that in A only the topology matters). The observed 
effect is quantum mechanical in nature: on the diagonal sites 
constructive interferences are particularly manifest. 

The phenomenon of fast CTQW transport through the net- 
work can also be observed when starting at any other node, 
where on short time scales we find a high probability of going 
to the "mirror" node. We define the "mirror" node of {jx,jy) 
to be (N + 1 — j x , N + 1 — j y ), so that the two nodes are 
related by inversion with respect to the center of the network. 

This selective behavior has to be contrasted with the one 
displayed by classical CTRWs. These describe namely sim- 
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pie diffusion, in which no exceptional transition probabilities 
between particular nodes occurs. 



may be contrasted to the classical CTRW, in which the limit- 
ing probability distribution is uniform, thus symmetric for all 
nodes. 



C. Limiting probabilities 

We continue by examining the situation at even longer 
times and focus in particular on the LPs Xk,j given by 
Eq.(|8j. Here, because we had already determined all eigen- 
values and eigenvectors, we took advantage of the structure 
of Eq.(|8}. From Eqs.0 and (|8}, and denoting the orthonor- 
malized eigenstates of the Hamiltonian by \q n ), such that 
E„l<Zn)(<zJ = 1. we find that 



^(fc|e- ;Ht | g „)(qJj) 



Yfc 7 = lim — dt 



= lim — I dt 



E e " <7Ant ( fe l9„)(g„b') 

n 

^2( k \Qn)(Qn\j)U\Q m )(Qm\ k ) 



I ^lim - / ^ e" l(A "- Am)7 * 



(17a) 



= ^2Sx n ,x m (k\q n )(q n \j)(J\q m )(q m \k). (17b) 

n,m 

We note that the integral in Eq.( ll7a> equals 1 if A„ = A m 
and otherwise, i.e., it equals S\ n ,\ m - Given that some eigen- 
values of H are degenerate, the sum in Ea.( ll7b> can contain 
terms belonging to different eigenstates \q n ) and \q m ). Equa- 
tion Mlb\ provides a numerically very efficient way of com- 
puting the Xk,j- Remarkably, we find that the \k,j depend 
in an unexpected way on the exact value of the size N of the 
finite network under study. Given these unexpected findings, 
which we report in the following, we cross-checked our evalu- 
ation method based on Ea.( ll7b> very carefully, by comparing 
it in selected cases to the direct evaluation of the integral in 
Eq.(|8}. In so doing, we fixed the upper integration limit to a 
very large value and verified that even larger values didn't lead 
to any changes in Xk,j- In a U cases we found that both numer- 
ical methods agree to very high precision. Thus, in general 
we prefered to work with Eq.l ll7bt , which is computationally 
much faster than Eq.(|8}. 

In Fig|8]we present Xk,m for CTQWs, in which the starting 
node is the middle node m. We display results for networks of 
sizes N = 9 and N = 15. Note that the Xk,m are symmetric 
about the initial middle node m, meaning that a node and its 
"mirror" node have the same limiting probabilities. This is in 
no way surprising, because this symmetry is already inherent 
in A and H. More remarkable are the patterns obtained. They 




Xk,m 
0.04 



FIG. 8: Limiting probabilities Xk,m to be at node k when starting at 
the middle node m for networks of sizes (a) N — 9 and (b) N = 15. 
The results are obtained from Eq. <17bt . 



Also when starting at a corner node c, we often find that 
the LPs for the starting node and its "mirror" node are equal. 
In Fig|9]we show the \k,c obtained by going from the corner 
node c = (1, 1) to the other nodes for networks of sizes N = 
5,N = U,N = 23, and N = 47. 

However, for some particular network sizes the distribu- 
tions of the LPs turn out to be asymmetric. For instance, for 
a network of size N — 15 the LP Xoc,c for the CTQW start- 
ing at node c to be at the opposite corner node oc is less than 
the LP Xc,c to be at the initial node. The same is true for 
the nodes along the edges of the network. Figure [lOl shows 
that such asymmetries occur for networks of the sizes N = 6, 
N = 15, N = 24, and N = 48 (the asymmetries are best 
seen by looking at Xc,c and Xoc,c)- Note that these asymme- 
tries occur for networks in which N is increased only by unity 
compared to networks which behave symmetrically, e.g., see 
Fig|9] The smallest network where we detected asymmetries 
in the distribution of the LPs has N = 6. The next ones we 
found for N = 12, 15, 18, 21, 24, 30, 36, ••• . 
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FIG. 1 1 : Differences between the LPs for CTQWs that start at c = 
(1, 1) to be at c, Xc,c, or to be at its "mirror" node oc = (N, N), 
Xocc, as a function of the network size N, for 1 < N < 60. 

The asymmetries are small and therefore not easy to spot in 
the global pictures displayed in Figs|5]and[TO] As illustrative 
examples, we choose prominent points in the network to show 
the asymmetries. An asymmetric LP distribution is particu- 
larly evident in the difference between Xc,c and Xoc,c- Thus, 
as an overview we present in Fig^2 as a function of A a plot 
of the (xc,c — Xoc,c)N 2 values obtained. Note that all A val- 
ues in Fig^2f° r which (xc,c — Xoc,c) 7^ are divisible by 3. 
However, the converse is not true, we find symmetric LP dis- 
tributions for the networks with A = 3,9, 27, 33, 39, • ■ ■ . The 
general A-dependence of Xoc,c{N) and Xm,m{N) is plot- 
ted in Fig^] In contrast to the classical limiting probability, 
which in all cases shows equipartition between all nodes, i.e. it 
is given by A , we find that Xoc,c 

(A) and Xm,m(N) decay 
nearly algebraically, namely as A~ 3 / 2 . 
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FIG. 12: iV-dependence of the LPs Xoc,c{N) and x™,m(N). For 
Xm,,m (N) only odd numbered networks are shown, since only those 
have a single central node m. 

Another striking feature of the CTQW is that the LPs dis- 
play quite regular patterns over the network. For the CTQW 
starting at the network's middle node m, the distributions of 
the LPs show a star-like pattern, see Fig|H] In all cases studied 
here, the LP distribution is such that its major fraction ap- 
pears to be distributed along lines diagonal and parallel to the 
network's edges and crossing each other at the initial node. 
Here again this implies that the transport is generated by con- 
structive and destructive interference. This might have conse- 
quences: On a regular, square network (or on a lattice for that 
matter) the application of the CTQW as a search algorithm 
is flawed, since its effectiveness is correlated with the initial 
site. That is, there is a high probability of finding a certain 



node which is "constructively" correlated with the initial one, 
but there is also a rather low probability of finding the oth- 
ers. Thus, although the topology of the square network has 
no exceptional sites, the transport through the lattice strongly 
depends on the initial condition. This relates directly to previ- 
ous studies, where it was found that quantum transport can be- 
come much slower than the classical one [8]. However, there 
it was also shown that if one starts in a superposition of states, 
the transport can get to be much quicker than in the classical 
case. A similar effect should also be observable here. For in- 
stance, if one starts in a uniform superposition of states along 
the baseline of the network, the CTQW can be mapped onto a 
one dimensional problem, as treated for instance in |7]. 

Furthermore, the fact that the LP distributions have espe- 
cially high peaks at the initial node and at its "mirror" node 
strongly recalls (in the spatially discrete version discussed 
here) the quantum mirage effects found in elliptic quantum 
corrals; thes e, again , can be related to wave interferences, see 
for instance 12711280 . A more detailed study of this effect will 
be published elsewhere 1 29] . 



IV. CONCLUSION 

We have studied numerically continuous-time quantum 
walks on finite networks topologically equivalent to square 
lattices. Furthermore, we compared our results to analytic ex- 
pressions obtained from the Bloch ansatz for networks with 
periodic boundary conditions. For these quantum walks, we 
have found that on short time scales a directed transport 
through the network takes place. In particular, when placing 
the initial excitation at one corner node, the walks propagate 
in a rather direct fashion along the diagonal to the opposite 
corner node. The transport is not ballistic, but is rather due to 
constructive quantum mechanical interferences. 

In the long time limit, we found that walks on networks 
of specific sizes A may show (in a totally unexpected way) 
asymmetric limiting probability distributions. This asymme- 
try manifests itself in the fact that the limiting probabilities for 
a CTQW to be at the initial node c = (1, 1) and at its "mirror" 
node oc = (A, A) differ. However, we were unable to find a 
way to predict which particular A values are related to such 
asymmetries. 

In general, the limiting probability distributions show pat- 
terns which depend on the starting site of the CTQW. This is a 
remarkable effect, which might have consequences for search 
algorithms based on CTQWs. Furthermore, we also found in 
all our calculations that the limiting probability distributions 
show strong peaks at the initial node and its "mirror" node. 
This effect resembles a discrete version of quantum mirages. 
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